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i The prime candidate of LIGO/VIRGO sources of gravitational waves is the spiral in of 

, black holes and neutron stars in compact binaries. While the early stages of the evolu- 

tion of compact binaries is computable from post-Newtonian calculations, prediction of 
their late stages requires large scale numerical simulation. A fully covariant and strictly 
hyperbolic formulation for numerical relativity is described, and illustrated in a one- 

O dimensional computation of a Gowdy-wave on the three-torus. This formulation allows 
foliations in full generality, in particular it poses no restriction on the lapse function. 
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1 Introduction 

> 

The advanced stages of construction of LIGO and VIRGO has stimulated some 
new developments in hyperbolic formulations for numerical relativitytxE3. These 
developments describe nonlinear gravitational wave motion in the late stages of 
spiral in of compact binaries of black holes and neutron stars. Binary black hole and 
neutron star systems remain today the prime candidates for astrophysical sources 
of gravitational radiation, and statistical calculations indicate their event rates to 
be. more than a few per year at the level of the sensitivity of the initial LIGO 
I , an. Other known burst sources with potentially detectable levels of gravitational 

^Jy radiation are the collapse of white dwarfs into neutron starsS, and secularly unstable 
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single neutron starsO. 



5— i ' 2 /(/) diagram for burst sources 
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The time evolution of the above mentioned three classes of burst sources of gravi- 
tational radiation can be illustrated (at the present level of our understanding) in 
an /(/)— diagram (Figure 1). Recall that compact binaries show an initial spiral in 
set by the Newtonian relationship 

(1) 

where / = d//dt. The frequency / refers to the base frequency of the gravitational 
wave-form, initially of the "chirp," and finally of the quasi-normal mode ring down. 
During the intermediate merger phase, / and / either follow by interpolation, or 
from / = ^<j> with <j> the orbital phase. Numerical relativity on black hole bina- 
ries specifically targets the merger phase, where post-Newtonian calculations break 
down. The present challenges are found to be the computation of stable, multiple 
orbits and the treatment of horizon boundary conditions. 
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GRAVITATIONAL WAVES FOR LIGO /VIRGO 
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FIGURE 1. Three cartoons of burst sources of gravitational waves in an /(/)-diagram. 
The top diagram shows the three phases of binary coalescence: the chirp (with the New- 
tonian result /(/) ~ / 11/3 ), the merger phase, and the ring down phase with quasi-normal 
mode frequency Jqnr ~ 865^^- Hz (see, e.g., Flanagan & Hughesu). The merger phase 
of two black holes and two neutron stars is expected to be distinct, but in as yet unknown 
ways. Sketched is the epmected trend with increasing spin S of the binary objects. Order 
of magnitude estimatesutllj of the radiated energy in the merger phase is also indicated. 
(A possible break up and disk formation in neutron star mergers has been by Zhugeu, and 
could also be included here.) The middle diagram sketches the expected gravitational 
wave form resulting from a dynamical instability in an (accretion induced) collapse sce- 
nario, such as the collapse of a white dwarf to a neutron star. The arrow indicates the 
direction of time in this trajectory. The first stage shows a rapid sweep of increasing fre- 
quency, followed by a phase of rapid damping due to gravitational radiation. This scenario 
assumes the ratio /3 of kinetic energy to gravitational energy to exceed the value which 
triggers & dynamical instability (/3 > 0.27). It has been suggested as a model for GRBs 
by Usovu. The lower diagram sketches a translation of results on the secular instability 
of single neutron stars, /3sec(= 0.1375) < j3 < f3([yn(~ 0-27), of Lai & ShapiroQ into the 
/(/)— diagram. The secular instability describes gravitational radiation driven transitions 
from unstable Jacobi (J u ) ellipsoids to either stable Mclaurin spheroids (M s ) or stable 
Dedekind ellipsoids {D a ). The Jacobi ellipsoids are rotation (fl) dominated (rel. vortic- 
ity), and the Dedekind ellipsoids are vorticity (ui) dominated (rel. rotation) stars. The 
trajectories are extended to the left, where the shape is in counter rotation relative to the 
fluid's vorticity. Note that the end point of the trajectory in the middle diagram connects 
to J u of the lower diagram. 
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3 Nonlinear wave equations for relativity 



Fully covariant, nonlinear wave equations have recently been obtained general rela- 
tivity in the tetrad approach. They are Yang-Mills type equations for the Utiyama- 
connections Lo a ^ v = (e M ) c V a (e ;y ) c (for a tetrad field {(e M ) a }^ = i) in the Lorentz 
gauge V a u> afll/ . In vacuo, they reduce to 

Ou afll/ - [w c , V a u> c ]^ v = 0, (2) 

were □ = V C V C is the £0(3, 1) gauged wave-operator. The full set of equations 
further involves the equations of structure, which can be integrated in the form of 

(e M ) a (T)=A/(T,ro)(e M ) c (ro)+ f A„"(t, s)d b N u (s)ds, (3) 

J T 

where the spatial coordinates have been suppressed. Here, the 

= (e M ) T (4) 

are the tetrad lapse functions, which are algebraically equivalent to the Hamil- 
tonian lapse and shift functions. Strict hyperbolicity is maintained with arbitrary 
N^, provided they define space- like foliations. There is full generality in this regard. 
By comparison, the otherwise closely related hyperbolic formulation of Abrahams 
et alr3 shows some restrictions on the Hamiltonian lapse function for maintain- 
ing strict hyperbolicity. (This difference was left unnoticed in the comparison by 
Fimfl). 

A numerical implementation is given in van Putten,E9 and applied to a nonlinear 
Gowdy-wave in a Gowdy T 3 -cosmology. Stable and accurate results are obtained 
with proper second-order convergence by comparison with analytical solutions. No 
indication of any "coordinate shocks" as reported by Alcubierre & Masso Ej have 
been observed in these numerical experiments. 
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